The isotope effect in the Hubbard model with local phonons 
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The isotope effect (IE) in the two-dimensional Hubbard model with Holstein phonons is studied 
using the dynamical cluster approximation with quantum Monte Carlo. At small electron-phonon 
(EP) coupling the IE is negligible. For larger EP coupling there is a large and positive IE on the 
superconducting temperature that decreases with increasing doping. A significant IE also appears 
in the low-energy density of states, kinetic energy and charge excitation spectrum. A negligible 
IE is found in the pseudogap and antiferromagnetic (AF) properties at small doping whereas the 
AF susceptibility at intermediate doping increases with decreasing phonon frequency wo- This IE 
stems from increased polaronic effects with decreasing luq. A larger IE at smaller doping occurs due 
to stronger polaronic effects determined by the interplay of the EP interaction with stronger AF 
correlations. The IE of the Hubbard-Holstein model exhibits many similarities with the IE measured 
in cuprate superconductors. 



I 

C 

O 

o 



> 

\o 
in 
q 

o 

oo 
o 



X 



Introduction. Isotope effect (IE) measurements [l[ 
played a significant role in establishing a phonon- 
mediated mechanism in conventional superconductors. 
One of the greatest successes of the Bardeen-Cooper- 
Schrieffer (BCS) theory of superconductivity was the 
ability to explain the IE experimental data. However, 
the IE in the high T c materials 0, 0, HI is unusual. At 
optimal doping the IE is small, leading many researchers 
to argue that phonons do not play an important role in 
the mechanism of high T c superconductivity and conse- 
quently that its origin must be purely electronic. Never- 
theless, the very large IE at small doping clearly shows 
that phonons play an important role in the physics of the 
cuprates. 

Some investigators have proposed the unusually large 
IE at small doping may be caused by the vicinity to a 
quantum critical point, defining a crossover from three- 
dimensional to two-dimensional (2D) physics as the dop- 
ing decreases |6|. Other studies claim that the IE in 
cuprates can be attributed to polaron formation Q- 
However, most theoretical investigations of the IE in 
models for high T c superconductors [7|, |8|, |f| start with 
approximations that assume a BCS-type mechanism for 
d-wave pairing and then study the effect of phonon char- 
acteristic energy on the superconducting properties. In 
contrast, the present study employs the dynamical clus- 
ter approximation (DCA) (lol 11 1 and treats, on equal 
footing, both the EP and electron-electron interactions 
in a 2D system in the form of the Hubbard-Holstein 
(HH) model. In previous work employing the DCA 
the single-band Hubbard model, a generally acce pted 



ccepted 



paradigm for the low energy physics of cuprates 
shows d-wave superconductivity driven by strong elec- 
tronic correlations. Furthermore, when combined with 
EP coupling, the superconducting temperature is signifi- 
cantly influenced by polaronic effects which are enhanced 
in the presence of antiferromagnetic (AF) correlations 
[H|,[l6| inherent in the underlying Hubbard model. Here, 
by polaronic effects we understand renormalization of the 



single-particle propagator, i.e. reduction of quasiparticle 
(QP) weight and charge carrier mobility due to EP cou- 
pling, and not necessarily the formation of small polarons 
or a narrow polaronic band. Thus far no investigation 
of the IE in systems with both strong electronic correla- 
tions and EP interaction has been conducted with cluster 
mean-field methods. 

For conventional superconductors the BCS theory pre- 
dicts that the isotope coefficient a = — 2d\J^ = 0-5, 
where T c is the superconducting temperature and luq the 
phonon frequency. The HH model with low frequency 
phonons, displays a positive and unusual IE. For small 
EP coupling the IE is small (a -C 0.5) while for val- 
ues of EP coupling large enough to drive the system in 
to a regime with strong polaronic effects a can be very 
large (a 3> 0.5) and a decreases with increasing dop- 
ing. The IE in the HH model is directly related with the 
enhancement of polaronic effects with decreasing ujq, as 
seen in the dependence of the single-particle quantities 
and charge susceptibility on changes in wo- A larger IE 
at smaller doping is due to increased polaronic effects in 
the presence of stronger AF correlations as described in a 
previous investigation The pseudogap and AF prop- 
erties display a negligible IE at small doping, whereas the 
AF susceptibility at intermediate doping increases with 
decreasing loq. The size, sign and doping dependence of 
the IE in the HH model is similar to experimental data 
in the cuprates 0, 17, HI- 

Formalism. The HH Hamiltonian reads 
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Here t is the nearest-neighbor hopping, U is the on-site 
Coulomb repulsion between electrons, ujq is the energy of 
the dispersionless optical phonons and g is the on-site EP 
coupling. The dimensionless EP coupling is defined as 
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A = 2g 2 / (2Mwp W) and represents the ratio of the single- 
electron lattice deformation energy E p — g 2 / (2Mujq) to 
half of the electronic bandwidth W/2 = At. The IE re- 
sults from the isotopic change of the ions atomic mass M 
and, since ljq oc M~? , A is kept fixed and loq varies. Re- 
sults are presented for small phonon frequency uo < 0.3t 
and for U — 8t, generally accepted or canonical param- 
eter regimes for studies of this type and in a range that 
may be appropriate for the cuprates [lij ]. 

The DCA [10|, a cluster mean-field theory that 
maps the original lattice model onto a periodic cluster 
of size N c — L 2 embedded in a self-consistent host, is 
used to study the HH Hamiltonian of fTJ). Correlations 
up to a range L c are treated explicitly, while those at 
longer length scales are described at the mean-field level. 
With increasing cluster size, the DCA systematically in- 
terpolates between the single-site mean-field result and 
the exact result, while remaining in the thermodynamic 
limit. Cluster mean-field calculations on the simple Hub- 
bard model [ll| or Hubbard model with local phonons 
rave been successful at capturing many key fea- 
tures associated with cuprate materials including a Mott 
gap and strong AF correlations, <i-wave superconductiv- 
ity and the pseudogap phenomenon. 

A quantum Monte Carlo (QMC) algorithm [20| modi- 
fied to perform the sum over both the discrete field used 
to decouple the Hubbard repulsion and the phonon field 
x serves as the cluster solver in the DCA. The present 
calculations are restricted to clusters of size N c — 4. Low 
temperature results for larger clusters are inaccessible 
due to the sign problem in the QMC calculation. The 
Maximum Entropy method [2l[ is employed to calculate 
the real frequency spectra. 

The superconducting properties are analyzed using the 
pairing matrix M = Txo- According to the Bcthc- 
Salpeter equation 



X = Xo + Xa^X = Xo(l - M)~ 



(2) 



the instability in the two-particle pairing propagator, 
X, indicating a possible transition to a superconduct- 
ing state, occurs when the leading eigenvalue of M 
reaches one. The corresponding eigenvector $ deter- 
mines the pairing symmetry. M is the product of 
the pairing interaction vertex T and the bare bubble 
Xo = GG. r can be seen as an effective paring in- 
teraction; xo is a convolution of the dressed single- 
particle propagators. The effect of EP coupling and luq 
on these two quantities can be studied by looking at 
Vd = Y,K,i^K>,^®d{KMnKM\K\iu>')$ d {K\i^ 
and P d = YsK^^diK'^XoiKiiuj). They are the re- 
spective projections of the interaction vertex T and xo 
on the subspace spanned by the (i-wave eigenvector $, as 
discussed in Ref [23| . 

Results. Previous DCA calculations show that at 
small doping large changes occur in the system properties 
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FIG. 1: (color online) T c versus uiq at 5% doping for A = 
0.3,0.5,0.6. Inset: T c versus A for uj = 0.3i, 0.2i, O.li. The 
transition to the polaronic regime (0.5 < A < 0.7) is accom- 
panied by an abrupt increase of the IE coefficient a. 



once A > A c « 0.5. The crossover (CRO) region defining 
the transition from large to small polaron regime, 0.5 % 
A ^$ 0.7 [22j is characterized by strong polaronic effects, 
a suppressed but finite T c and a large IE effect. Outside 
the CRO region, for small EP coupling a negligible IE 
effect is found, whereas for A > 0.7, T c vanishes due 
to formation of small and heavy polarons. The IE for 
different values of A is shown in Fig. [TJ When A = 0.3, 
T c is nearly independent of ujq, and just slightly smaller 
than T c of the Hubbard model without phonons (see the 
inset). At A = 0.5 we find a significant reduction of T c 
with decreasing phonon frequency c^o although the IE 
coefficient a w 0.13 is much smaller than the BCS value 
of 0.5. However, for A = 0.6, i.e. inside the CRO regime, 
the IE is very large once uiq < 0.3t. Here the IE coefficient 
isasj 1.41. 

We find that the IE is decreasing with increasing dop- 
ing. Fig. [5] (a) shows the logarithm of T c versus the log- 
arithm of t^o at different dopings when A = 0.5. To cal- 
culate the IE coefficient we use a linear fitting of ln(T c ) 
versus 2 ln(w ) for u < 0.4i, the slope giving the esti- 
mate of a for this study. The IE coefficient dependence 
on doping is shown in Fig. [5] (b). The decrease of a 
with increasing doping is in agreement with the notion 
of significant interplay between the AF correlations, EP 
interaction and superconductivity [HI]. Given stronger 
AF correlations at smaller doping, the polaronic effects 
become more pronounced and hence the IE is larger. It 
would be interesting to study the doping dependence of 
the IE for larger values of A (e.g. A = 0.6). Unfortu- 
nately at large A and intermediate doping (i.e. « 15%) 
the charge susceptibility of the system becomes extremely 
large implying strong charge fluctuations that make it 
very difficult to stabilize the calculation at low tempera- 
ture. 
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FIG. 2: (color online) (a) T c versus oj at 5%, 10%, 15% , 20% 
and 25% doping for A = 0.5, on a logarithmic scale. A linear 
fitting of ln(T c ) versus 21n(o;o) provides an estimation of IE 
coefficient a from the slope, (b) a versus the filling n at 
A == 0.5. The error bars are of the order of the symbol size. 
a decreases with increasing doping. 
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FIG. 3: (color online) (a) The d-wave pairing interaction Vd 
and (b) the d-wave projected bubble Pd versus temperature 
for different u>o, at A = 0.6 and 5% doping. While the pairing 
interaction increases with decreasing cjo, the decrease of Pd 
results in an overall suppression of T c with decreasing u)o- 



In regard to superconductivity, the EP coupling gives 
rise to two competing effects [16|, an increase of the ef- 
fective pairing interaction and a renormalization of the 
single-particle propagator. The first favors superconduc- 
tivity while the latter opposes it. Keeping the EP cou- 
pling fixed and decreasing loq results in a similar com- 
petition; i.e., an increase of the effective pairing interac- 
tion and a stronger renormalization of the single-particle 
propagator leading to the same overall effect, a reduction 
of T c . As discussed in Ref. 23j, the d- wave projected 
quantities Vd and Pd provide, respectively, a measure of 
the effective d-wave pairing interaction and renormaliza- 
tion of the single-particle propagator relevant for <i-wave 
superconductivity. At low temperature as ojq decreases 
Vd increases while Pd decreases, as shown in Fig. |3l 

In the CRO regime the IE on the superconducting T c is 
accompanied by a significant IE in the low energy DOS 
and in the kinetic energy of the system. For A = 0.6 
when wo decreases, the DOS in an energy window of or- 
der w loo around the chemical potential (chosen zero in 
FigH] (a)) is suppressed and the kinetic energy at low 
temperature increases with decreasing u>o, as shown in 
Fig. 2] (b). Like the projected bubble Pd discussed in 
the previous paragraph, both the DOS and the kinetic 
energy describe single-particle propreties of the system 
and are strongly linked to effects on charge carrier mo- 
bility and effective mass. Their dependence on uj shows 
that the polaronic effects become stronger as the phonon 
frequency is reduced, a fact understood by noting that 
decreased luq moves the system closer to the adiabatic 
regime where the CRO region shrinks to A c and charge 
carriers should eventually self-trap when A > A c [3]. 

The magnetic properties are also influenced by the 



phonon energy. At small doping and at half filling the 
magnetic susceptibilities show very small luq dependence, 
on the order of the error bar, even in the CRO regime (not 
shown). However, the IE on the Neel temperature TV 
can be significant near the doping range where TV — > 0. 
For example, as can be seen in Fig. [4] (c) at 15% dop- 
ing, the AF susceptibility increases as u>o decreases. The 
pseudogap temperature T* measured as the characteris- 
tic temperature for suppression of low energy spin excita- 
tions is virtually independent of uiq. It is interesting that 
when T* is measured as the temperature at which the 
pseudogap first begins to appear in the DOS, the same 
conclusion can be drawn. Moreover, one can see in Fig. [4] 
(a) that the magnitude of the pseudogap measured by the 
peak-to-peak distance increases weakly as ujo decreases. 

The charge excitation spectrum shows a strong depen- 
dence on ojq. In the CRO regime the local (i.e. r = 0) 
charge excitation spectrum displays a narrow peak in 
agreement with the assumption that the system tends 
toward formation of a polaronic band [HI]. Moreover, 
and also in agreement with that assumption, the width 
of the peak is proportional to u>o as shown in Fig. Q] (d). 
The low-energy weight of the charge spectrum increases 
as luq decreases and explains the stronger fluctuations 
in the electronic density and the difficulty in fixing the 
doping level as ujq decreases at low temperatures. 

Discussions. While the HH Hamiltonian shows an IE 
in many respects similar to the one seen in cuprates, as we 
will discuss in the next paragraph, a direct comparison 
with experiments requires more realistic models which 
consider the symmetry of the relevant phonon modes. 
Since the relevant phonons are the ones associated with 
oxygen vibrations, the explicit inclusion of the oxygen 
orbitals in the model might be important as well. 

However, we find many similarities between the IE in 
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FIG. 4: (color) The black (red, blue) lines represent wq = 
0.3t (wo = 0.2t, luq = O.lf) results for A = 0.6 while the 
green ones show data for A = 0. (a) DOS iV(w) and (b) 
kinetic energy Ekin = y~] k e(k)nk, a at 5% doping. The low 
energy DOS and the kinetic energy gain are suppressed as wo 
decreases, c) The AF susceptibility x sp at 15% doping. The 
AF correlations are enhanced when wo decreases, d) The local 
dynamic charge susceptibility Xch( r = 0, w)/w at 5% doping. 
The width of low energy peak in the charge excitation spectra 
is proportional to ojq. 



cuprates and our results in the CRO regime. The exper- 
iments in cuprates report a decrease of a with increasing 
doping [3J , a negligible IE in the pseudogap and AF prop- 



181 and a large IE in the spin 
25] indicative of an increase of 



erties at small doping 17, 
glass freezing temperature 
AF correlations at finite doping with decreasing phonon 
frequency. The large IE seen in the penetration depth 
shows strong dependence of the Cooper pairs effective 
mass on phonon frequency 0, 0, [26|. Although we do 
not address the superconducting state, this may be un- 
derstood with our findings of strong loq dependence of the 
charge carriers effective mass. The competition between 
the enhancement of the pairing interaction and the renor- 
malization of the single-particle propagator (discussed in 
Fig. [3]) may be relevant for explaining the negligible IE 
in the superconducting T c accompanied by a significant 
IE in the charge carrier effective mass seen in the opti- 
mally doped cuprates @, since the former is a result of the 
competition and the latter is a consequence of only the 
single-particle renormalization. The experiments find an 
extremely large IE around 8% doping reportedly linked 
with stripe formation 0,H3|. Stripes cannot be addressed 
with the small 2x2 clusters since the the stripe phase or- 
der parameter is incommensurate with the cluster. How- 
ever, the strong enhancement of charge fluctuations due 
to EP coupling and the dependence of charge susceptibil- 
ity on loq suggest a large IE on stripe formation if larger 
clusters were considered. 

Due to the nature of cluster-mean field approximation 



which neglects spatial correlations outside the cluster 
range, in the underdoped region where the phase fluc- 
tuations (PF) are strong, the polaronic effects' influence 
on the PF is underestimated. This leads to the expecta- 
tion of an even larger IE at small doping and a stronger 
doping dependence of a in calculations which manage to 
capture longer range correlations. This can be under- 
stood by noting that the PF coherence temperature is 
given by the where n s is the superconducting fluid 
density and m* is the effective mass of the pairs [28j |. 
Thus increased effective mass due to the decrease in ojq 
amplifies the PF and reduces T c . 

Conclusions. This letter presents results of calcula- 
tions of the two-dimensional Hubbard model with Hol- 
stein phonons using DCA and QMC specifically address- 
ing the IE in light of recent experiments on high T c mate- 
rials, including the cuprates. At small EP coupling there 
is a negligible IE, but for larger EP values, in the CRO 
regime, a large and positive IE effect is found on the su- 
perconducting T c that becomes stronger at small doping. 
A significant IE is seen also in single-particle quantities, 
like the low energy DOS and the kinetic energy, as well as 
the charge excitation spectrum. This indicates that the 
IE on T c is directly related with the enhancement of ef- 
fects such as suppression of QP weight and charge carrier 
mobility with decreasing ujq. A larger IE at smaller dop- 
ing is a consequence of an increase in these same effects 
in the presence of strong AF correlations. A negligible IE 
is seen in the pseudogap and AF properties at small dop- 
ing whereas the AF susceptibility at intermediate doping 
increases with decreasing phonon frequency u>q. 
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